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The upper bound on the mass of the light CP-even Higgs boson of the minimal 
supcrsymmetric model (MSSM) depends on the supcrsymmetric particle spectrum 
via radiative loop effects. At present, complete one-loop results and partial two- 
loop results are known. Some simple analytic approximations have been obtained 
which are quite accurate over a large portion of the MSSM parameter space. Based 
on these results, I examine how accurately one can predict the upper bound on the 
mass of the lightest MSSM Higgs boson. 



1 Introduction 

Low-energy supersymmetry [0 provides the most compelling framework for 
electroweak physics, in which the electroweak symmetry breaking is generated 
via the dynamics of an elementary scalar Higgs sector. The scalar boson masses 
are kept light (of order the electroweak symmetry breaking scale) due to an 
approximate supersymmetry in nature. Supersymmetry is broken at the TeV 
scale or below, and this information is transmitted to the scalar sector, thereby 
generating electroweak symmetry breaking dynamics at the proper scale. 

The simplest model of low-energy supersymmetry is the minimal super- 
symmetric extension of the Standard Model (MSSM). In this model, the Higgs 
sector consists of eight degrees of freedom made up from two complex weak 
scalar doublets of hypercharge ±1 respectively 0. Supersymmetry requires 
that the hypercharge —1 [+1] Higgs doublets couple exclusively to down-type 
[up-type] fermions, respectively. After minimizing the Higgs potential, the 
neutral components of the Higgs doublets acquire vacuum expectation values 
(vevs) with (Hi) — Uj/v2. The model possesses five physical Higgs bosons: 
two CP-even scalars, h° and H° (with m^o < mjjo), a CP-odd Higgs scalar A 
and a charged Higgs pair H ± . As usual, I define tan (3 = v^/vx and normalize 
v 2 = v\ + v\ = im^/g 2 — (246 GeV) 2 . Due to the form of the Higgs-fermion 
interaction, the third generation quark masses are given by mi, — h^vxj \[2 and 
m t = h t V2/V2, where h q (q — t, b) are the corresponding Yukawa couplings. 

The tree-level physical Higgs spectrum is easily computed Its most 
noteworthy feature is the upper bound on the light CP-even Higgs scalar: 
rriyfi < mz\cos2[3\ < niz- The maximum tree-level upper bound of mz is 
saturated when one of the vevs vanishes (and m^o > mz)- It is convenient to 
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consider a limiting case of the MSSM Higgs sector where v\ = 0. For finite 
this limit corresponds to m& = 0, which is a reasonable approximation.^ In the 
vi = model, the Higgs sector degenerates to a one-doublet model with: 

V mggs = m 2 &<S> + i A ($t$) 2 , X ^\{g 2 +g' 2 ). (1) 

The supersymmetric constraint on the value of A is a consequence of the fact 
that the MSSM Higgs quartic couplings originate from the D-term contribu- 
tions to the scalar potential. The squared-mass of the light CP-even Higgs 
boson of the v± = model is given by m 2 = Xv 2 = m 2 z . 

2 Upper bound of m h o in the MSSM 

The upper bound of m^o < mz will be modified by radiative corrections . 
In order to obtain the radiatively corrected upper bound of m^o, it suffices 
to compute radiative corrections in the V\ = model. Let us focus on the 
real part of the neutral scalar component: $° = (v + h)/y/2. The bare Higgs 
potential takes the following form: Vniggs = toh + ^(m^oh 2 + 0(h 3 ), where 

to = v + ml] , (m| )o = ™o + | A o«o . ( 2 ) 

and the subscript indicates bare parameters. We also introduce 

("4)o = i(5 2 + .9o>o- (3) 

The on-shell renormalization scheme is defined such that mz and m^o are phys- 
ical masses corresponding to zeros of the corresponding inverse propagators. 
Let the sum of all one-loop (and higher) Feynman graphs contributing to the 
Z and h° two-point functions be denoted by iAzz(q 2 )g fJ ' L ' + iBzz(q 2 )q fl q L ' and 
—iAh^q 2 ), respectively, where q is the four- momentum of one of the external 
legs. The physical masses are given by: 

m z = Ol)o + Rc A zz (m%), (4) 
m 2 = (m 2 o)o + Re A hh (m 2 h o) . (5) 

Since v is the vev of the scalar field at the true minimum of the potential, we 
require that the sum of all tadpoles must vanish. That is, 

t + A h (0) = 0, (6) 

"In practice, it is sufficient to take Vi -C V2, and then fix the value of /if, to be consistent 
with the observed fe-quark mass. 
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where — iAh(0) is the sum of all one-loop (and higher) Feynman graphs con- 
tributing to the h° one-point function. Combining eqs. (|^)-(^), we obtain 

m 2 h o = mi + Re [A hh {m%) - A zz (m%)] [ Ao - ±(. 9o 2 + 3o 2 )] v 2 . (7) 

This result is accurate at one- loop order, since we have put m^o = rriz and 
vq = v on the right hand side where possible. 

Naively, one might argue that eq. (Q) can be simplified by using the su- 
persymmetric condition Ao = jig 2 . + g' 2 )- However, this is correct only if a 
regularization scheme that preserves supersymmetry is employed. Of course, 
the physical quantity m 2 must be independent of scheme. Consider two dif- 
ferent regularization schemes: dimensional regularization (DREG) and dimen- 
sional reduction Q (DRED). Renormalized couplings defined via (modified) 
minimal subtraction in these two schemes are called MS and DR couplings, 
respectively. DREG does not preserve supersymmetry because the number of 
gauge and gaugino degrees of freedom does not match in n =/= 4 dimensions. In 
contrast, DRED preserves supersymmetry (at least at one and two-loop order). 
In DRED [DREG], bare quantities will be denoted with the subscript D [G]. 
Then, the supersymmetric condition holds in DRED: 

\D-l{9l+9%) = 0. (8) 

We now demonstrate that the above relation is violated in DREG. First, 
the gauge couplings of the two schemes are related as follows Q 

gl = 9% + ^, g'3 = g'3- (9) 

For the Higgs self-coupling A, the relation between the two schemes is derived 
by considering the one-loop effective potential (in the Landau gauge), V = 
y(o) _|_ where is the tree-level scalar potential and V' 1 ' is given by: 



=-— ^StrX 4 ($) 



if -In 



M 2 {>S>) 



M 2 



(10) 



In eq. (|l0|), K is a scheme-dependent constant (see below), A4 2 (§) denotes 
the squared-mass matrix as a function of the scalar Higgs fields (i.e., the cor- 
responding tree-level squared-mass matrices are obtained when <!> is replaced 
by its vev), and the divergences that arise in the computation of the one-loop 
integrals in 4 — 2e dimensions appear in the factor A = 1/e— 7^ + ln(47r) [where 
7e is the Euler constant] . We have also employed the notation 

Str {• ■ •} = J2 C t (2J t + f )(-l)" +1 {• • ■}« , (11) 
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where the sum is taken over the corresponding mass matrix eigenvalues, in- 
cluding a factor Cj which counts internal degrees of freedom (e.g., charge and 
color) for all particles of spin Ji that couple to the Higgs bosons. 

In DRED, K = 3/2, independent of particle i in the sum [eq. (0)]. The 
fact that particles of different spin yield the same constant K is an indication 
that DRED preserves supersymmetry at one-loop. In DREG, K = 3/2 for 
spin and spin-1/2 particles, while K — 5/6 for spin-1 particles. However, the 
effective potential (expressed in terms of bare parameters) must be independent 
of scheme. Comparing the DREG and DRED computations, it follows that 

^dv 4 - ^ (|) (6mk + 3m|) = ±A G « 4 - (§) (6m 4 w + 3m 4 ) , (12) 



which yields 

Ad 



ff 4 (m| + 2m 4 ,) 

Ag "I fi/i 2 4 ■ \ 16 > 



Combining the results of eqs. (||) and ( |l3|) gives the DREG result 



m + g'S) = -tt^-^ K + §"4) • (14) 



Thus, in computing the one-loop corrected Higgs mass [eq. (^)] in DREG 
[DRED], one must use the relation between Ao and j(go+g'o) given be eq. (14) 
[eq. (H)]. One can check that this difference is precisely compensated by the 
difference in DREG and DRED that arises in the computation of the vector 
boson loop contributions to the one-point and two-point functions. Henceforth, 
we shall always use the DRED scheme, in which case [|| 

m 2 h o = m| + Re [A hh (m 2 z ) - A zz (m 2 z )] - . (15) 

v 

Although the loop functions above are individually divergent, all divergences 
precisely cancel in the sum and yield a well-defined one-loop result for m^o. 

The method described above [resulting in eq. ([if])] is sometimes called 
the diagrammatic method [Sin since one explicitly evaluates the one-point 
and two-point functions by standard Feynman diagram techniques. A second 
method for computing m h o, called the effective potential technique, is often 
employed in the literature This is not an alternate "scheme", but 

simply another way of organizing the calculation. Consider the DRED one- 
loop effective potential introduced above (with K = 3/2). In the sum V = 
y(o) _j_ the DR scheme consists of absorbing the factor of A into the 

bare parameters (mo, Ao and $o), which converts them into DR parameters. 
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Renormalized quantities (such as the effective potential or the n-point Green 
functions) will be denoted with tildes in the following. These are computed in 
the Landau gauge; the divergent piece A is removed by DR subtraction and 
the bare parameters are replaced by renormalized DR parameters. Finally, the 
DR parameters are related to physical parameters. 

We proceed as follows. First, we minimize the renormalized effective po- 
tential by setting the first derivative equal to zero. This condition yields: 

[\\v 2 + m 2 ] v + A h {0) = 0. (16) 

In eq. (|l6|), the first term on the left hand side arises at tree-level, while the 
second term is a consequence of the the fact that the nth derivative of 
evaluated at the potential minimum, is equal to the scalar n-point function 
evaluated at zero external momentum |Tcj . The second derivative of the effec- 
tive potential, denoted by (m^ )eff, is similarly given by: 

{m 2 h o) e ff = m 2 + %\v 2 + A hh (0) . (17) 

We may use the DRED relation [eq. (||)] , which is also satisfied by the renor- 
malized DR parameters, to eliminate A. The DR Z-mass parameter is given 
by (m|)j5R = j(g 2 + g' 2 )v 2 - Combining the above results yields: 

Ko)eff = (m|) M + A hh (0)-^1. (18) 

In the literature, (ni^o)eff is sometimes used as the approximation to the 
one-loop-improved Higgs squared-mass. However, this is not a physical pa- 
rameter, since it depends on an arbitrary scale that is introduced in the DR 
subtraction scheme. To obtain an expression for the physical mass, which cor- 
responds to the zero of inverse propagator, we note that (rn^ )eff has been 
computed using the two-point function evaluated at zero external momentum. 
Thus, the physical Higgs squared-mass is given by: 

mlo = (m|o) e ff + Re A hh (m 2 h0 ) ~ A hh {0) . (19) 

Likewise, we must convert from (rnz)-^ to the physical Z squared-mass. This 
is accomplished using a result analogous to that of eq. (^) , which guarantees 
that mz corresponds to the zero of the inverse Z propagator: 
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^z = (m 2 z ) m + ReAzz(m 2 z). (20) 
Combining eqs. (|l8|)-(p0|), we end up with 

m\ a = m| + Re [l^(m|) - A ZZ {m 2 z )\ - . (21) 

Not surprisingly, we have reproduced the diagrammatic result [eq 
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3 Leading Logarithms and Renormalization Group Improvement 



When the loop functions in eq. (^) are computed, one finds that the most sig- 
nificant contributions grow logarithmically with the top squark masses. (Terms 
that are logarithmically sensitive to other supersymmetric particle masses also 
exist.) Over a large range of supersymmetric parameter space, the radiatively 
corrected Higgs mass can be well approximated by just a few terms. On the 
other hand, if the logarithms become too large, then the validity of the pertur- 
bation theory becomes suspect. However, in this case the leading logarithms 
can be resummed using renormalization group (RG) techniques 012) 



We begin with a one-loop analysis. Consider an effective field theory ap- 
proach Jl2[ , and assume for simplicity that supersymmetry breaking is char- 
acterized by one mass scale, Msusy, which is assumed to be large compared 
with mz- At scales /i < Msusy, the Higgs potential takes the form: 



V=\m 2 {n)[h{n)Y + \\{n)[h{^ 



(22) 



Letting h — > h + v with m 2 (jj,) < 0, the Higgs mass is given by 



m\ {n) = \(fi)v 2 (n) 



(23) 



Since the effective theory is supersymmetric only for /x > Msusy, we impose 
the supersymmetric boundary condition [see eq. dq)] : 



A(Msusy) = \ [g 2 (Msusy) + 5' 2 (M SU sy)] 



(24) 



Scale dependent parameters satisfy renormalization group equations (RGEs). 
For \i < Msusy, the Standard Model RGEs are relevant: 



ft 



lv 



dX 



i 



d In fi 2 
dlnv 2 



16tt 2 
1 



din /i 2 167T 2 
d(g 2 +g' 2 ) 



6A 2 



W + (g 2 +g' 2 ) 2 ]-2j2N cf h 



din fi 2 



96tt 2 



(8iV g -43) 3 4 + (f7V 9 + lV' 



2A 7 „, 



(25) 



where hj = \/2mf/v, N g = 3 is the number of fermion generations, and 
N c f — 3 [1] when / runs over quark [lepton] indices. 

It is instructive to solve the RGEs iteratively to one-loop, by ignoring the 
/j dependence on the right hand sides in eq. (25). Incorporating the boundary 
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condition [eq. (|24|)], the solution for X(mz) is given by]] 
A(m z ) = i(<? 2 + 5 ' 2 )(M SUSY )-/3 A lnf^% SY 



(.9 2 + 5' 2 )(-z) + (|/W - A) In (Mk^X.) 

\ m z / 



1 f„2 
4 



(26) 



Finally, using eq. (|23|), we identify the physical Higgs mass by evaluating 
mf i0 (fi) at /i = mz and taking v(mz) = 246 GeV. We know from the previous 
section that this is not strictly correct. However, at the one-loop leading 
logarithmic level, this procedure is accurate, and we end up with: 

K„) 1LL = m% + (i/? ff2+9 , 2 - f3 x )v 2 In (MkfA } (27) 

\ m z J 

where the subscript ILL indicates that the result is only accurate to one-loop 
leading logarithmic order. To obtain the full one-loop leading logarithmic 
expression, simply insert the results of eq. (25) into eq. ( p7j ) [in f3\ one can 



consistently set A = \{g 2 + g' 2 )]- We have checked p|,|l2f that the above result 
matches precisely with the diagrammatic computation [eq. Jl5|) l in the limit 
of Msusy 3> where Msusy characterizes the scale of supersymmetric 

particle masses (taken to be roughly degenerate). The dominant term at one- 
loop is proportional to mf and arises from the piece of /3 a proportional to hf . 
Inserting (3\ = —3hj/8ir 2 with h t = y/2m t /v into eq. (p7|), one obtains j] 

(4)ilt = m| + -g^ In ( *^L) . (28) 



' w 



\ rn? t J 



The subscript 1LT indicates that this is the leading mf piece of (m 2 i0 )iLL- 
However, the additional terms in (m 2 t0 )iLL are numerically significant as we 
shall show at the end of this section. 

Thus, we see that given the RG functions, no additional diagrammatic 
computations are needed to extract the full one-loop leading logarithmic con- 
tribution to the Higgs mass. Thus the RG-approach provides a useful short 
cut for extracting the leading one-loop contributions to the Higgs mass. Of 
course, if the leading logarithms are large, then they should be resummed to all 



6 Onc subtlety consists of the proper way to run down from trit to mz, since below 
fi = mt, the electroweak symmetry is broken. I will ignore this subtlety here although it can 



: addressed; see ref. 



The lower scale of the logarithm in this case is mj* (and not rn l z ) since this term arises 
from the incomplete cancelation of the top quark and top squark loops. 
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orders. This is accomplished by computing the RG- improvement of the exact 
one-loop result as follows. Let (to? )irg = A(mz)« 2 (raz), where X(mz) is 
obtained by numerically solving the one-loop RGEs. Write the exact one-loop 



result as: 



(m? ) 



ILL 



Iinl, where (m^ )iNL is the result obtained 



by subtracting the one-loop leading logarithmic contribution from the exact 
one-loop result. Clearly, this piece contains no term that grows logarithmi- 
cally with Msusy- Then the complete one-loop RG-improved result is given 
by m 2 h0 = (m^o)iRG + ("^o)inl- 

The RG technique can be extended to two loops as follows jl3|. For sim- 
plicity, we focus on the leading corrections, which depend on at = h 2 /4ir and 



f/47T, 



we work in the approximation of = g = g' = and 



A <C h t . (All two- loop results quoted in this section are based on this approx- 
imation.) The dependence on the strong coupling constant is a new feature 
of the two-loop analysis. We now solve the one-loop RGEs by iterating twice 
to two loops. In the second iteration, we need the RGE for h 2 , which in the 
above approximation is given by 

fa -jA- hi = [Ik* - 8 5 2 ] hi . (29) 



dln/x 2 167r 2 

This iteration produces the two-loop leading double logarithm 0], and yields 



AN = W+9' 2 )- 



3hf(m t ) 



In 



M 2 

ivi SUSY 



A- 
hi 



in 



^SUSY A 



8tt 2 \ mi ) \' u hf J \ m 2 J 

(30) 

Next, we must incorporate the sub-dominant two loop effects. Only three 
modifications of our one-loop analysis are required (in the limit of = g = 
g' = and A <C ht). First, we need only the ht and g s dependent parts of the 
two loop contribution to (3\. That is, (3\ is modified as follows pi 



1 



[30^-32^ 



4 „2 " 



(31) 



(167T 2 ) 2 

Including this into the iterative solution of the RGEs adds a two-loop singly 
logarithmic term to the result of eq. (p0|). Second, we must distinguish between 
the Higgs pole mass (denoted by m^o with no argument) and the running Higgs 
mass evaluated at m t . Using the results of Sirlin and Zucchini [jlEj , 



m 2 h0 



Am 2 w \{m t ) 



1 



1 



(2± 

V 7T 



(32) 



Third, we make use of the relation between v (mt) and v = Am^/g 



v 2 (m t ) 



(?) 



(33) 



Using the above results, we end up with 
3 5 2 



2 2 



87T 2 TO^ 



m t (m t ) In 



M 2 

JW SUSY 



Iv 



Phi 



hi 



/" ^SUSY 



V 7T 



3 /a* 

8 V 7T 



(34) 



where /i t = ht(m t ) and m t (m t ) = h t (m t ) v(m t )/V2. Numerically, the two- 
loop singly logarithmic piece of eq. ( |34| ) contributes about 3% relative to the 
one-loop leading logarithmic contribution. 

Let us compare this result with the two-loop diagrammatic computation of 
ref. |l6| ]. In order to make this comparison, we must express eq. (|34|) in terms 
of the top quark pole mass, mt. The relation between m* and the running 
top-quark mass is given by |l7| , |l8t 



mt = m t {m t ) 



3 V 7T / 2 V 7T 



(35) 



where m t (m t ) is the MS running top-quark mass evaluated at m t .|] Inserting 
the above result into eq. (Bjj) yields: 



3 g 2 mf 



In 



( -^SUSY 



4a, 



7o 



13 



In 



M 2 

JW SUSY 



(36) 



This result matches precisely the one obtained in ref. |lq| in the limit of 
Msusy S> mz- Note that the numerical contribution of the two- loop singly- 
logarithmic contribution in eq. (|3^ ) is about 10% of the corresponding one-loop 
contribution. Clearly, the use of the running top quark mass [as in eq. ([34]) ] 
results in a slightly better behaved perturbation expansion. 

Finally, we can employ a very useful trick to make our results above even 
more compact. The two-loop doubly-logarithmic contribution can be absorbed 
into the one-loop leading-logarithmic contribution by an appropriate choice of 
scale for the running top-quark mass. Specifically, using the iterative one-loop 
leading-logarithmic solution to the RGEs for h t and v yields 



1 



3a t 
16tt 



hi 



rn. 



(37) 



d We caution the reader that ref. [ha] defines mt(mt) = ht(mt)v / y/2, which differs slightly 
from the definition of mt {mt ) used here. 
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If we choose the scale [it = VmtMsvsY to evaluate the running top-quark 
mass in eq. (p4[) , we end up with: 



m 2 h0 



3ff 2 



57T^TO 



■mt(Pt) In 



M 2 

■'"SUSY 

mf(jUi) 



1 + lfM_iW 
3 V tt / 16 V tt / 



(38) 



One can check that the sum of the terms in the brackets deviates from one 
by less than 1%. Thus, in practice, the two-loop singly-logarithmic contri- 
bution can now be neglected since it is numerically insignificant. That is, 
one can incorporate the leading two-loop contributions by simply inserting 
the running top-quark mass evaluated at [i t = y/rrit Msijsy into the one-loop 
leading- logarithmic expression for m^ • 
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Figure 1: The upper bound to the mass of the light CP-even Higgs boson of the MSSM 
is plotted as a function of the common supcrsymmetric mass MgusY ( m the absence of 
squark mixing). The one- loop leading logarithmic result [dashed line] is compared with the 
RG-improved result, which was obtained by a numerical computation [solid line] and by the 
simple recipe described in the text [dot-dashed line]. Also shown are the leading m* result 
of eq. (teSf) [higher dotted line], and its RG-improvement [lower dotted line]. The running 
top quark mass used in our numerical computations is mt(mt) = 166.5 GeV. 

Fig. |l| illustrates the results of this section. We display the results for m^o 
based on five different expressions for the light CP-even Higgs mass. Case (i) 
corresponds to the one- loop leading m\ result, (m^ )iLT [eq. (J28|)] . In case (ii) 
we exhibit the full one-loop leading logarithmic expression, (to^T) ill [eq. (|27|)]. 
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In case (iii), we consider (to? )irg obtained by solving the one-loop RGEs 
numerically. Finally, case (iv) corresponds to the simple recipe proposed above, 
in which we evaluate (m^ )iLL by setting m t to the running top quark mass at 
the scale fit- For completeness, we also include case (v), where we apply the 
same recipe to (to^o)ilt- 

The following general features are noteworthy. First, we observe that 
over the region of M S usy shown, {m h o) 1RG ~ {m h o) lhL (m t (fi t ))- Second, the 
difference between (m^o)iLL and (m^o)iRG is non- negligible for even moderate 
values of -Msusy; neglecting RG- improvement can lead to an overestimate of 
m h a which can be as large as 10 GeV (for Msusy > 2 TeV, the deviation 
grows even larger). Finally, note that although the simplest approximation, 
(m^o)iLT, reflects the dominant source of radiative corrections, it yields the 
largest overestimate of the light Higgs boson mass. 

4 Additional Complications: Supersymmetric Thresholds 

In the analysis of the previous section, we assumed that all supersymmetric 
particle masses were roughly equal and substantially larger than mz- To ac- 
count for a non-degenerate supersymmetric spectrum, we must recompute the 
RGEs in steps starting from fi = Msusy and ending at mz- Every time the 
threshold of a supersymmetric particle is passed, we integrate it out of the 
theory, and determine a new set of RGEs for the new effective theory. Eventu- 
ally, when we pass below the lightest supersymmetric threshold, we regain the 
RGEs of the Standard Model given in cq. (H|). We can solve iteratively for 
A(mz) as we did in the previous section, but now using the more complicated 
set of RGEs. Explicit formulae can be found in refs. jl2| and |l3| |. 

However, the above procedure fails to incorporate the effects of squark 
mixing. Since the most important contribution to the Higgs mass radiative 
corrections arises from the incomplete cancelation of the top quark and top 
squark loops, it is important to examine this sector more closely. First, we 
define our notation. The physical top squark squared-masses (in the v\ = 
model) are eigenvalues of the following two 2x2 matrix 

/ M^ + m 2 -m 2 z (±- e t sin 2 W ) m t A t 
\ m t A t Mfj + m 2 — m 2 z e t sin 2 9w 

where e t = 2/3 and Mq, My, A t are soft-supersymmetry-breaking parameters. 

We shall treat the squark mixing perturbatively, assuming that the off- 
diagonal squark squared-masses are small compared to the diagonal terms. | 

e Formally, we assume that (Mf - M$)/(M% + Mf) < 1, where M\, Mf are the top 
squark squared-masses. Thus, we demand that mtAt/M^ VSY <C 1. 
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The perturbative effect of squark mixing is to modify the supersymmetric re- 
lation between the Higgs quartic coupling and the gauge couplings [eq. (p4|)]. 
Such modifications arise from one loop squark corrections to the Higgs quartic 
self-coupling via: (i) corrections to the scalar two-point function on the exter- 
nal legs; (ii) triangle graphs involving two trilinear Higgs-squark-squark inter- 
actions and one quartic Higgs-Higgs-squark-squark interaction; and (iii) box 
graphs involving four trilinear Higgs-squark-squark interactions |19|. Then, 
eq. ( pi| ) is modified to: 

(40) 



A(Af SUSY ) = \(g 2 + g' 2 ) + SX 2 + 5X 3 + SX 4 , 
where the 5\ arise from the three sources quoted above. Explicitly, 



5\ 2 = 
SX 3 = 
5\ 4 = 



-3(5 2 



1^-2 ■ 



32tt 2 



u) 



, ;2 2 v AhiA 2 t h{M 2 Ql M 2 
^- 2 hiAtg(M 2 ,M 2 



AihiB{M 2 Q ,M v 

g 2 h 2 A 2 Pt (M 2 ,M 2 ) 



(41) 



where 



B(a,b) 

h(a, b) 

f(a,b) 

g(a,b) 
Pt(a,b) 



1 



(a - b) 2 
1 



\ (a + b) 



In - 



(a-b) 
1 

(a - b) 2 
f(a,b) + 



-b 



In 



a b 



In 



\b 









:?)] 


a — o 





2et sin 2 9\y(a — b)g(a, b) . 



(42) 



For simplicity, consider the case of Mq — Mjj = Msusy- Using B(a, a) = 
l/6a, h(a, a) = 1/a, /(a, a) = — l/2a and g(a, a) = — l/6a 2 , eq. (J40|) becomes: 



A(M SU s Y ) = i(5 2 + .9' 2 ) + 



3h}A 



8tt 2 M| usy 



A 2 t 



1 2 M 2 



(43) 



Note that the correction term due to squark mixing has a maximum when 
At — \/6-Msusy- This relation is often called the maximal mixing condition, 
since it corresponds to the point at which the one-loop radiative corrections to 
m^o are maximal. 
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Using the new boundary condition, we may repeat the analysis of the 
previous section and recompute mf iQ . At one loop, the effect of the squark 
mixing is simply additive. That is, the modification of m^ due to squark 
mixing at one loop is given by: (Am^imix = (6X2 + <5A 3 + SX^v 2 . At two- 
loops, we solve for X(mz) by iterating the RGE for A(/i) twice as in the previous 
section. However, the boundary condition for A(Msusy) bas been altered, and 
this modifies the computation. The end result is 



, . 2N 3g 2 mfAf ( Af 

(Am 2 h ) mix = o » IJ 1- 



4? A , (M, 



r 2 

"? . 
(44) 

i.e., (Am?)niix acquires a logarithmically-enhanced piece at two loops. In this 
approximation, the maximum in (Am^) m i X at A t = VQMsvsy is not shifted. 
However, this method does not pick up any non-logarithmically-enhanced two- 
loop terms proportional to A t . To obtain such terms, one would have to 
perform a two-loop computation in order to find the necessary two-loop terms 
that modify the boundary condition [eq. (f40|)]. 

It is again possible to absorb the two-loop singly-logarithmic term into the 
one-loop contribution, (Am^)i m i X , by an appropriate choice of scale for the 
top-quark mass. The end result is quite simple: 



_ 3g 2 m t(M SUSY )A 2 ( A 2 
(Am ' J — 8^M| USY V 12Mf us J' 



That is (Am^) m i X = (Am^)i m i x (mt(/Lt~)), where the appropriate choice of scale 
in this case is /i~ = Msusy- The difference from the previous case [where 
/it — V^^Msusy] arises due to the extra factor of 2 multiplying the two-loop 
singly- logarithmic term in eq. (Q) [compare this with eq. (|30|)]. Physically, 
fi~ = AfsuSY corresponds to the scale at which the squarks decouple and the 
boundary condition [eq. (f43| )l is modified due to squark mixing. 

To illustrate the above results, we compare in Fig. |^ the value of m^o as 
a function of A t based on the five cases exhibited in Fig. 0. Specifically, the 
effects of (Am^) m j X are included at the one-loop level in cases (i) and (ii) , 
while cases (iv) and (v) make use of the improved result given by eq. ([45|). 
In the full RG-improved result [case (iii)], the RGE for A(/i) is computed 
numerically using the modified boundary condition [cq. (ff(i|)1. We see that 
(mlo)iRG — ("4o)iLL(m t (/i t )) + (Am^ )i mix (m t (/I t )). Thus, once again a sim- 
ple recipe provides an excellent approximation to the numerically-integrated 
RG-improved result over the entire region of the graph. Note that the maxi- 
mal value of myfi occurs for \A t \ ~ 2.4Msusy- The solid or dash-dotted line 
provides our best mass estimate, and we conclude that m/jo < 125 GeV for 
-^susy < 1 TeV. Similar results were also obtained by Carena et al. M. 



13 



a 



150 




RG-improved (full) ( ) 

RG— improved (approx) ( 

one— loop leading log ( 


-) 


140 




.leading m\ results ( ) 




130 


_ / ~~~~ 

/---■. 
J_ 
{ 


X 

\ 

\ 

\ / 


, " " x - 

/ 






• . \ ' . / 
■ \ '■ . / 


120 


// 


■ . N /. 
=x ■ \ /. ' 




110 




1 TeV ^ "rnf le = 


175 GeV 



-3-2-10 1 2 3 

A t /M SUSY 



Figure 2: The upper bound to the mass of the light CP-even Higgs boson of the MSSM 
plotted as a function of j4t/MsuSY- Squark-mixing effects are incorporated as described in 
the text. See the caption to Fig. ll 



During the past year, two groups have computed the A t dependence of m^o 
at the two-loop level. Ref. pOfl has performed a diagrammatic two-loop compu- 
tation which includes all terms of 0(a s ), as a function of tan (3. Ref. uses an 
effective potential approach to extend the computation of ref. |l6[ and compute 
directly the two-loop squark mixing contributions in the v± = model. These 
results show that the A t dependence of m h o is modified slightly at two loops: 
the maximal squark mixing point occurs at A t ~ 2A/susy, a value somewhat 
below the result noted above. Moreover, the value of m h o at maximal squark 
mixing is slightly higher than the one shown in Fig. ^ for MgusY = 1 TeV, 
the maximal value of m^o is found to be close to m h o ~ 130 GeV. Presumably, 
these results are due to genuine two-loop non-logarithmically enhanced terms 
proportional to a power of A? /M| TTSV . An important check of the calcula- 
tions presented in refs. j^o) and [|lj would be to explicitly verify the two-loop 
logarithmically-enhanced contribution exhibited in eq. ([l4|). 

5 Conclusions 

I have described in detail the theoretical basis for the computation of the upper 
bound of the mass of the light CP-even Higgs boson of the MSSM. It suffices 
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to consider the limiting case of v \ = which considerably simplifies the anal- 
ysis. I explained how one can use renormalization group methods to provide 
a short-cut for obtaining the leading one-loop and two-loop contributions to 
myfi . These methods can also be generalized to the full MSSM Higgs sector at 
arbitrary tan (3. Further details and references can be found in ref . Jl^| . 

As a result of the work by many groups during this past decade, we be- 
lieve that the predicted value of m h o as a function of the MSSM parameters 
is accurately predicted within an uncertainty of a few GeV. Simple analytic 
formulae provide an excellent representation of the known results over a large 
range of the MSSM parameter space The present partially known two- 

loop information is essential to this conclusion and provides confidence that 
there are no surprises lurking in some corner of the supersymmetric parameter 
space. Some clarification is still needed to understand more completely the 
dependence on the squark mixing parameters. 
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